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On the Mass of the neutron star in Cyg X-2 
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ABSTRACT 

We present new high resolution spectroscopy of the low mass X-ray binary Cyg X- 
2 which enables us to refine the orbital solution and rotational broadening of the 
donor star. In contrast with Elebert et al (2009) we find a good agreement with 
results reported in Casares et al. (1998). We measure P = 9.84450 ± 0.00019 day, 
K2 = 86.5 ± 1.2 km s _1 and V smi — 33.7 ± 0.9 km s _1 . These values imply q = 
M 2 /M 1 = 0.34 ± 0.02 and Nh = 1.71 ± 0.21 M (for i = 62.5 ± 4°). Therefore, the 
neutron star in Cyg X-2 can be more massive than canonical. We also find no evidence 
for irradiation effects in our radial velocity curve which could explain the discrepancy 
between Elebert et al's and our K2 values. 
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1 INTRODUCTION 

Dynamical studies in low mass X-ray binaries (LMXBs here- 
after) offer a promissing route to test the equation of state of 
nuclear mat ter. Soft equations of s tate are not stable above 
1.6 Mq (e.g. iBrown fc Bcthc 1994) and hence finding a neu- 
tron star more massive than this limit would be a major ad- 
vance in our knowledge of nuclear matter physics. LMXBs 
and, in particular, accreting millisecond pulsars (AMPs) are 
expected to harbour the most massive neutron stars because 
of the sustained accretion of matt er during their long lifes 
(|van den Heuvel fc B itzaraki 1995). Unfortunately, dynam- 
ical studies are hampered by (1) the overwhelming accre- 
tion luminosity in persistent LMXBs (which swamps the 
donor star's spectrum) and (2) the extreme faintness of 
the companion star in transient AMPs during quiescence 
jD'Avanzo et al.ll2009r ). 

A promis sing route to overc o me th ese limitations was 
proposed by ISteeghs fc Casaresl l|2002i '). X-ray irradiated 
donor stars can be betrayed by the detection of high ex- 
citation fluorescence NIII and CIII lines within the Bowen 
blend at ~ 4634-50 A. The Doppler shift of these narrow lines 
traces the orbit of the donor star and, therefore, can yield 
dynamical constraints in persistent LMXBs and transient 
AMPs during outburst. As a result of this strategy some 
good candidates for massive neutron stars hav e been found, 
namely X1822-371 jMunoz-darias et al.1 120051 ') and Aql X-l 
|Cornelisse et al.ll2007n ~ 

In a few long period (P > 1 day) systems, the donor 
is an evolved star which can be spectroscopically detected 
over the irradiated accretion disc. Cyg X-2 represents one of 



thi s exceptional cases. I ts first orbital solution was presented 
bv lCowlev et~ai] i|l979h who show that the donor star is an 
A5-F2 III orbiting the neutron star in 9.843 days with a pro- 
jected velocity K2 = 87 ± 3 km s~ . Al most 20 years later, 
this was refined bv lCasares et all (|l998l ) (C98) who find an 
A9 donor with orbital parameters P olb = 9.8444 ± 0.0003 d 
and K2 = 88.0 ± 1.4 km s~ . This work also reports the first 
determination of the rotational broadening of the donor's 
absorption features (V sin i = 34.2 ± 2.5 km s" 1 ) which, in 
turn, implies a binary mass ratio q = M2/M1 = 0.34 ± 0.04. 
This value is remarkable as it implies a peculiar donor star, 
very undermassive for its spectral type an d the probable 
outcome of an intermediate mass binary llKing fc R itter 
1 19991 : iPodsiadlowski fc Rappaportl l2000l : iKolb et al.ll2000l ). 
But the implications for the compact object's mass are also 
remarkable. The orbital solution, combined with inclination 
constraints i — 62.5 ± 4° derived through ellipsoidal fits 
to UBV l ight curves gives a neut ron star mass of 1.78 ± 
0.23 M Q |o rosz fc Kuulkerslll999l ) and, hence, a good can- 
didate for a massive neutron star. This r esult has been re- 
cently challenged by lElebert et ail (|2009l l (E09) who mea- 
sure K 2 = 79 ± 3 km s" 1 and hence Mi=1.5 ± 0.3 M©. 



Here we present new high resolution spectroscopic ob- 
servations of Cyg X-2 with the main aim of refining the ro- 
tational broadening of the donor star and update the orbital 
parameters. In contrast with E09, we find a good agreement 
with previous results reported in C98 which give support to 
the presence of a massive neutron star in Cyg X-2. 
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2 OBSERVATION AND DATA REDUCTION 

Cyg X-2 was observed on the nights of 25-26 July 1999 us- 
ing the Utrech Echelle Spectrograph (UES) attached to the 
4.2 m William Herschel Telescope (WHT) at the Observa- 
torio del Roque de Los Muchachos. Ten 1800-3600s spec- 
tra were obtained with the E31 echelle grating and 2Kx2K 
SITel detector, covering the spectral range AA5300-9000 A. 
We selected a 1" slit and a factor 2 binning in the spectral 
direction, resulting in 10 km s _1 resolution. ThAr arc images 
were observed each night for the purpose of wavelength cal- 
ibration. The final wavelength calibration was verified using 
the 01 skylines at A5577.34 and A6300.304 and variable ve- 
locity offsets were found between 0.2-5 km s _1 . These were 
corrected from every individual spectrum. 

Eleven additional spectra were obtained on the nights 
of 31 July 1999 and 9 July 2000 with the ISIS double-arm 
spectrograph on the WHT. Here we employed the 1200B 
grating on the blue arm and different central wavelengths 
resulting in wavelength coverages within AA3550-6665. A 1" 
slit was selected yielding spectral resolutions in the range 
34-54 km s _1 . Frequent comparison CuAr+CuNe arc lamp 
images were taken every night for wavelength calibration. 
This was tested against sky lines and it was found to be 
accurate to within 4 km s _1 . These small offsets were never- 
theless corrected from the individual spectra. The ISIS red 
arm was always centered redwards of 7600 A but is not 
used in this paper due to its lower spectral resolution and 
the paucity of absorption lines, which result in larger errors 
in the cross-correlation analysis. The red arm spectra were 
mainly taken for the sake of abundance analysis and will be 
reported elsewhere (Gonzalez Hernandez et al. 2009). A log 
of the observations is presented in Table [T] 

For the purpose of radial velocity analysis we observed 
the stellar template HR 114 (A7 III) using all different in- 
strumental configurations. This star has an intrinsic V sin i 
of 21 km s _1 i.e. comparable to the rotational broadening 
of the donor star in Cyg X-2 (C98). Therefore, to refine our 
previous rotational broadening we decided also to observe 
the F3 V star HR 6189 with UES, which has a reported 
upper limit V sini < 15 km s _1 . 

All the images were processed following standard de- 
biasing and flat-fielding, and the spectra subsequently ex- 
tracted using conventional optimal extraction techniques in 
order to opt imize the signal-to-noise ratio of the output 
l|Hornelll986h . 



3 REVISITING THE ROTATIONAL 
BROADENING 

In C98 we measured the rotational broadening of the donor 
star's features in Cyg X-2 using 25 km s _1 resolution spectra 
and found Usini = 34.2 ± 2.5 km s _1 . Here we revisit this 
determination using the 10 km s _1 UES spectra. We broad- 
ened the template star HR 6189 fro m 5 to 50 km s _1 in steps 
of 1 km s _1 , using a Gray profile (|Gravlll992r ) and a limb- 
darkening coefficient e = 0.5 appropriate for our spectral 
type and wavelength range. The broadened templates were 
rectified to the continuum (usig a low-order spline fit) and 
multiplied by factors / < 1, to account for their fractional 
contribution to the total light. These were subsequently sub- 
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Figure 1. A section of the Doppler corrected average of the Cyg 
X-2 UES spectra in the rest frame of the donor star (top), to- 
gether with the F3V template HR 6189 broadened by 33.7 km 
s — 1 (bottom). 

tracted from the Doppler corrected average of Cyg X-2, ob- 
tained using the orbital solution given in Sec. 4. The Cyg X-2 
average was also rectified through fitting a low order spline 
to the continuum, after masking out the main emission and 
atmospheric/IS absorption lines. The optimal broadening, 
based on a \ 2 test on the- residuals, is found for 34.6 ±0.1 
km s" 1 . A potential source of systematics is the assumption 
of continuum limb-darkening coefficient in the computation 
of the rotational profile. Absorption lines in late-type stars 
are expected to have smaller core limb-darke ning coefficients 
than the continuum (|Collins fc Truaxll 19951 ) and, therefore, 
assuming the continuum limb-darkening coefficient could 
bias the result. This was tested by repeating the same anal- 
ysis using zero limb-darkening as a conservative lower limit 
and obtain V sin i = 32.8 ±0.1 kms" 1 . Therefore, we decide 
to adopt the mean of the two limb-darkening values as a safe 
estimate of the true rotational broadening in Cyg X-2 i.e. 
V sin i = 33.7 ±0.9 kms" 1 . This is in excellent agreement 
with our determination in C98. 



4 UPDATED ORBITAL SOLUTION AND 
MASSES 

We rectified the 21 individual spectra by subtracting a low- 
order spline fit to the continuum, after masking out the 
main emission and absorption features. The UES spectra 
were subsequently rebinned into a uniform velocity scale of 

5 km s" 1 pix -1 and the ISIS spectra into 27 km s~ pix - . 
Radial velocities were extracted through cross-correlation of 
every individual spectrum of Cyg X-2 with its corresponding 
A7III template HR 114, observed with identical instrumen- 
tal setup. The template spectra were previously broadened 
to 33 km s _1 to match the width of the donor photospheric 
lines (see previous section) . Cross-correlation was performed 
in the spectral regions free from emission and telluric ab- 
sorption features. For consistency, new radial velocities were 
extracted from the old database (C98) using a contempora- 
neous observation of the same template HR 114, obtained for 
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the purpose of spectral classificatio n. Radial velocit i es wer e 
extracted following the method of iTonrv fc Davis! l| 19791) , 
where parabolic fits were performed to the peak of the cross- 
correlation functions, and the uncertainties are purely sta- 
tistical. The final errorbars were multiplied by a factor 2 so 
that the minimum reduced \ 2 °f a sinewave fit model is 1.0. 
This yields the following parameters 



7 = -209.9 ±1.4 kms -1 



P = 9.84450 ± 0.00019 d 



To = 2451387.148 ± 0.018 



K% = 86.5 ± 1.2 kms" 1 

where To corresponds to the inferior conjunction of the 
donor star. All quoted errors are 68 per cent confidence. The 
systemic velocity 7 has been corrected from the radi al veloc- 
ity o f HR 114, that we take as -10.2 ±0.9 km s" 1 |Wilson I 
1 1953! ). Fig.[2]displays the radial velocity points folded on the 
orbital period together with the best sine fit solution. New 
radial velocities are marked with open circles whereas solid 
circles indicate velocities from the old 1993-1997 database 
(C98). Our K2 velocity is consistent with C98 within 1-a 
but not with E09 who finds K 2 = 79 ±3 kms -1 . To ilustrate 
this point we also plot in dashed-line style the best sinewave 
fit for K2 = 79 km s -1 . This has xt — 3 for 60 degrees of 
freedom and hence it is far less significant than our K2 value 
ijLampton et al.lll976h . 

E09 have suggested that the difference between the two 
K2 values could be caused by different levels of X-ray irra- 
diation between the two observing epochs. The detection of 
Hell A4686 emission from the companion by E09 certainly 
indicates that the star is irradiated. However, it remains to 
be seen whether irradiation is sufficient, not only to pump 
chromospheric Hell emission, but also to modify the surface 
distribution of the photospheric absorption lines. In prin- 
ciple, irradiation can quench absorption lines from the in- 
ner hemisphere of the donor star, leadin g to an increase in 
the observed K2 l|Wade fe Hornelll988l ). It might be pos- 
sible that, by a chance coincidence, the E09 database was 
obtained during an episode of lower X-ray activity than our 
data and we have looked for this using RXTE/ASM (All 
Sky Monitor) data. Contemporaneous X-ray observations 
are available for the second half of the C98 campaign (since 
1996, when RXTE was launched) our new data from 1999- 
2000 and the entire E09 database but comparable levels of 
X-ray activity are found, with the X-ray flux Fx oscillating 
between 35 and 50 ASM counts s~ . Only one velocity point 
in C98, obtained on the night of 5 Aug 1996, was taken dur- 
ing a dip in the X-ray light curve of ~18 counts s _1 . Interest- 
ingly, its orbital phase is <f> — 0.63, almost identical to that 
of another point obtained on 3 Aug 1997 (4> = 0.67), when 
Fx — 50 counts s" 1 i.e. almost a factor 3 higher. Despite 
the difference in X-ray flux both velocities are consistent 
with our best orbital solution within 1— a. And one should 
note that phase 0.65 is close to an orbital quadrature, when 
velo city distortions from a circular orbit should be largest 
(e.g. iDavev fc Sniithll 19921 ). This strongly suggests that the 



effects of X-ray irradiation in the radial velocity curve are 
unimportant. 

As a matter of fact, irradiation will distort the radial 
velocity curve from a simple sine wave, introducing a ficti- 
tious eccentricity which should be measurable. Therefore, we 
have also attempted to fit e lliptical orbits to our database, 
following iFriend et al.l (|l990l ). Our best fit yields null eccen- 
tricity (e = 0.004 ± 0.019) and a larger reduced \ 2 than a 
simple circular solution, another indicati on that irradiation 
is neg ligible at this level. Furthermore, lOrosz fc Kuulkersl 
|l999l ) find no evidence for excess light at phase 0.5 in their 
optical light curves nor C98 observe any significant change 
of spectral type with orbital phase. These two results also 
suggest that X-ray irradiation is not enough to explain the 
discrepant K2 values. The lack of irradiation signatures is 
probably due to the fact that the donor star is ho t and the 
orbital separation large lOrosz fc KuulkeTsl J 19991 ). The X- 
ray flux received by every surface element on the companion 
star is too small to produce any disturbance in the radial ve- 
locity curve or light curve, despite the near Eddington X-ray 
luminosity of Cyg X-2. 

Looking at the radial velocity curve of E09 (shown in 
fig. 4) we note a large scatter in the velocity points near 
the phase 0.25 quadrature. Two datapoints seem to lie sys- 
tematically lower than the rest by ~20 km s _1 and this is 
certaintly dragging the K2 velocity to lower values. The au- 
thors admit that only one arc spectrum was obtained for 
most of the nights. Although sky lines were used to correct 
for instrumental offsets, the fact that the strongest sky line 
OIII A5577 lies at the edge of their spectral range may in- 
troduce some systematics in the offset correction. Therefore, 
we believe that the value reported by E09 might be affected 
by problems with the wavelength calibration which is a criti- 
cal issue given their low spectral resolution of ~160 km s . 
In any case, new high resolution observations obtained at 
the two orbital quadratures are clearly required to further 
constrain K2 and confirm our result. 

Since the donor star is filling its Roche lobe and syn- 
chronized, we can combine our updated K2 and V sin i values 
to constrain the binary mass ratio through 



Vsini = K 2 (1 +q) 



0.49 q 2/3 



(1) 



0.6 g 2 /s+ln(i + g i/3) 

ijHorne. Wade fc Szkodvlll986h which leads to q = 0.34 ± 
0.02. The revised mass function is thus /(Af) = 
Misin 3 i/(1 + qf = PKI/2-kG = 0.66 ± 0.03 and hence 
Mi sin 3 i = 1.19 ± 0.0 6 M . Assuming i = 62.5 ± 4° 
|Orosz fc Kuulkerslll999l ) we find Mi = 1.71±0.21 M . The 
error budget is clearly dominated by the uncertainty in the 
inclination angle and thus ellipsoidal model fits to new light 
curves are urgently needed to better constraint the stellar 
masses. 



5 SUMMARY 

We have revisited the determination of the system parame- 
ters in Cyg X-2 with 21 new high-resolution spectra obtained 
during 1999 and 2000. The new solution does not support 
the conclusions of E09 who claim for a significantly lower 
value for the radial velocity semiamplitude of the donor 
star. Instead, our results confirm previous determinations 
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Figure 2. Radial velocity curve folded on the ephemeris listed 
in Sect. 4 with the best sine-wave solution overplotted. New UES 
and ISIS velocities from 1999-2000 are marked with open circles 
whereas solid circles indicate velocities from the 1993-1997 ISIS 
database. The dashed line depicts the orbital solution of E09. 



reported in C98. The discrepancy cannot be explained by 
X-ray irradiation because the sinusoidal shape of the radial 
velocity curve is not disturbed. In particular, we find (i) no 
evidence for orbital eccentricity and (ii) no significant de- 
viations between two velocity points at phase ~ 0.65, when 
X-ray flux varies by a factor ~3. Our refined orbital parame- 
ters are P = 9.84450 ±0.00019 days, K 2 = 86.5 ±1.2 km s -1 
and V sin i = 33.7±0.9 km s~\ which lead to q = 0.34±0.02, 
Mi sin 3 i = 1.19 ± 0.06 M fl . Assuming i = 62.5 ± 4° from 
lOrosz fc Kuulkersl (ll999T ) leads to Mi = 1.71±0.21 M© and, 
therefore, the possibility that Cyg X-2 harbours a neutron 
star more massive than canonical. 
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Table 1. Log of the observations. 



Date 


Instrument 


Wav. Range 


Exp. time 


Dispersion 


Resolution 






AA 


(s) 


(A pix- 1 ) 


(km s _1 ) 


25-26/07/1999 


UES 


5300-9000 


3200, 3x(1800, 2700, 3600) 


0.11 


10 


31/07/1999 


ISIS 


5885-6665 


3x1800 


0.45 


34 


9/07/2000 


ISIS 


3545-4460 


6x1800 


0.45 


44 


9/07/2000 


ISIS 


4350-5250 


100, 1800 


0.45 


54 



